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RESUME. On demontre que le rapport signal a bruit, si important en theorie de I'information, devient 
sans objet pour des communications numeriques ou la demodulation s'effectue selon des techniques 
nouvelles d'estimation rapide. Calcul operationnel, algebre differentielle, algebre non commutative et 
analyse non standard sont les principaux outils mathematiques. 

ABSTRACT. The signal to noise ratio, which plays such an important role in information theory, is 
shown to become pointless for digital communications where the demodulation is achieved via new 
fast estimation techniques. Operational calculus, differential algebra, noncommutative algebra and 
nonstandard analysis are the main mathematical tools. 
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symboles, modulation, demodulation, bruits, estimation, rapport signal a bruit, calcul operationnel, 
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1 . L'auteur dedie ce travail a Claude LOBRY, chercheur aussi brillant qu'original, et ami fidele, qui lui 
a beaucoup appris, I'analyse non standard par exemple. 
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Extended english abstract 

Introduction 

The symbol to be transmitted is modulating a carrier, which is assumed to be a solution 
of a linear differential equation with polynomial coefficients. Most signals utilized in 
practice, Uke J2,]„kc sin(w,,t + i^,,), A,,, uj^, (pi G M, sine t = do satisfy this 

property. New algebraic estimation techniques ET\ l23l permit to achieve demodulation 
even with very "strong" corrupting additive noises and therefore to question the impor- 
tance of the signal to noise ratio, which is playing such a crucial role in information theory 
(see Il43ll44l and l2l|3]|4]|7m0l). Operational calculus, differential algebra and nonstan- 
dard analysis are the main mathematical tools. 

Identiflability 

Let ko{Q) be the field generated by a finite set Q — {9i, . . . , 6g} of unknown parameters, 
where ko is a ground field of characteristic 0. We are utilizing the classical notations of 
operational calculus ||48]| . Introduce the differential field |8, 41 1 fc(s) of rational functions 
in the indeterminate s over the algebraic closure k of ko{Q), with derivation Any 
signal X, X ^ 0, is assumed to satisfy a homogeneous linear differential equation with 
coefficients in k{s) and therefore to belong to a Picard-Vessiot extension |8, 41] of k{s). 
Write k{s) [j^] the noncommutative ring of linear differential operators with coefficients 
in k{s). The left fc(s)[J^] -module spanned by x and 1 is a fc(s)-vector space of finite 
dimension n + 1, n > 0. It yields the minimal non necessarily homogeneous linear dif- 
ferential equation Q where the polynomials p,qo, . . . ,qn in k[s] are coprime. Introduce 
the square matrix 9Jl of order N + M + 1, where the line, < ^ < ^ + is ©. 
Techniques stemming from Wronskian determinants ||8] |4T1 demonstrate that the rank of 
9Jl is + M. It follows that the coefficients of the polynomials p,qo, . . . ,qn are pro- 
jectively linearly identifiable 112111231 . Consider as a particular case x = where the 
polynomials p,q G k[s] are coprime. Then the coefficients of p and q are also projectively 
linearly identifiable. 

Perturbations and estimators 

Assume that the unknown parameters Q are linearly identifiable 112111231 . With and addi- 
tive perturbation w, we obtain the estimator ^ where 21, and 53, £ are respectvely g y. g 
and g X 1 matrices, such that the entries of 21 and 05 belong to spanj,^j^^j_d_j (1, x) and 

those of £ to span^[_d.j (w), where R is the localized ring I29II fco(0)[s](fco[s])^^- More- 
over det(2l) 7^ 0. It is always possible to obtain an estimator which is strictly polynomial 
with respect to -j, i.e., where the rational functions in s are polynomials in without con- 
stant terms. As in [11] it yields in the time domain, if x is an analytic function. Formula 
dTji where c is a constant, [0, t] is the estimation time window, the divisor 5{t) is an ana- 
lytic function such that (5(0) = 0, [6*1] e(i) is the estimated value of 9^ at time t. 

Noises 

We are considering two types of perturbations, which are noises in the sense of IfTTII : 
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- The first noise, which is zero-mean, is a finite sum X^nmic sin(f7ii + (pi) where the 
frequencies fti > are unlimited. 

- Let *N, *R be the nonstandard extensions [|42J of N, R. Replace [0, 1] C M by the 
hyperfinite set ||42l I = {0, . . . , 1}, where TV e *N is unhmited. A zero-mean 
white noise is a function w : I ^ *M, l > w{l) = An(i), where 

- A G *M is a constant, such that is limited, 

- the 71 (i) are independent zero-mean random variables with a normalized co variance 

1. 

The estimator (|7| yields "good" values for the unknown parameters for any limited values 
of the amplitudes Ai, A, and even for some unlimited values of them. 
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1. Introduction 

Le rapport signal a brui^ que Ton retrouve dans les formules de la theorie de 1' in- 
formation, telle qu'elle s'est imposee depuis Shannon (voir f??'/??) et, par exemple, dans 
la vaste litterature sur le sujet, |l2][3j|4]|7]|40|), est un ingredient fondamental pour definir 
la qualite des communications. Le but de ce travaijl est de demontrer qu'une nouvelle 
approche de 1' estimation rapide et du bruit (voir 1 11 1 et sa bibliographie) rend ce rapport 
sans objet dans un certain cadre numerique. Revoyons, done, le « paradigme de Shan- 
non ». Le symbole a transmettre (voir, par exemple, Il28ll40l ) module une porteuse z{t), 
solution d'une equation differentielle lineaire a coefficients polynomiaux : 

fmie 

La plupart des signaux utilises en pratique, comme une somme trigonometrique finie 
At sin(ti;ti + (yS, ), un sinus cardinal ou un cosinus sureleve , A^, cj^, 

tpi, uj G M, verifient une telle equation, qui se traduit dans le domaine operationnel (cf. 
IBHI ), pour t > 0, par 

^a.(-^Ki = J(s) (1) 

finie 

oil / e C[s] est un polynome dont les coefficients dependent des conditions initiales en 
< = 0. La demodulation revient, alors, a estimer certains des coefficients de ([Hi. On y 
parvient, ici, grace a des techniques algebriques recentes (cf. Il2ni23]| ). 

Un bruit, selon [11], est une fluctuation rapide, que Ton definit de fa9on efficace et ele- 
gante grace a 1' analyse non standarc^. Les calculs du §|4]sont effectues avec une somme 
finie de sinusoides a tres hautes frequences et un bruit blanc, dont la definition non stan- 
dard clarifie F approche usuelle des manuels de traitement du signal. lis demontrent la 
possibilite d'obtenir de « bonnes » estimations avec des bruits « tres forts », c'est-a-dire 
de « grandes » puissances, fait confirme par des simulations numeriques et des expriences 
de laboratoire (voir |l20l |30l |33l |3l [38l |45l |46 , 47 1). Les imperfections, inevitables en 
pratique, proviennent de I'implantation numerique des calculs, notamment de celui des 
integrates (voir ll30l[33l ). des interferences entre symboles (voir Il2l l28l[39ll40l et leurs bi- 
bliographies), et du fait que les bruits ne sont pas necessairement centres (voir a ce propos 
le § 3.2.2 de [TT!). 

Calcul operationnel et algebre differentielle aux §|2]et[3] analyse non standard au §|4] 
sont les principaux outils mathematiques. 



1. Rappelons au lecteur peu au fait que toute transmission physique de signal est perturbee, 
notamment par du « bruit ». L'extraction des informations utiles malgre ces alterations est un but 
essentiel en traitement du signal et en theorie des communications. Que Ton pense par exemple aux 
codes correcteurs d'erreurs, ou, comme les mathematiciens le savent, la theorie des codes en blocs 
est d'une grande richesse algebrique. 

2. Voir [13] pour une version preliminaire. 

3. Voir aussi [32] . 
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Remarque 1 Ajoutons pour le lecteur etranger, voire hostile, a 1' analyse non standard^ 
qu'il est loisible de la remplacer par des considerations « classiques ». On y perdrait en 
concision et, a notre avis, en intuition. 

Remarque 2 Avec des signaux analytiques par morceaux (le sens du mot analytique est 
celui de la theorie des fonctions et non pas, ici, celui usuel en traitement du signal ( cf. 
^ 13911401/ )). qui ne satisfont pas d' equations differentielles connues a I'avance, on utilise, 
selon les memes principes algebriques, des derivateurs numeriques afenetres glissantes 
pour obtenir les estimations ( voir A3 5]/ et M7\ 1791 12 7^. leurs exemples et leurs bibliogra- 
phies). On ne peut, alors, esperer les memes resultats que precedemment. 

Remarque 3 La possibilite de liens entre theorie de I' information et mecanique quantique 
a ete examinee par divers auteurs ( voir, par exemple, ^ |5] \24]l ). Rappelons a ce propos 
que I'approche du bruit en UIV a deja conduit a une tentative nouvelle de formalisation 
du quantique Iil2^ . qui sera completee grace a un resultat remarquable et tout recent, du 
a Charreton 

Remarque 4 Les techniques d' estimation evoquees plus haut ont permis des avancees 
notables en automatique^ lineaire ou non. 

Remarque 5 Des travaux en cours portant sur I 'ingenieriefinanciere devraient egalement 
demontrer I' applicabilite de notre point de vue a cette disciplin^ 



Remerciements. L'auteur exprime sa reconnaissace a O. Gibaru (Lille), M. Mboup (Pa- 
ris) et a tous les membres du projet INRIA-ALIEN pour des echanges fructueux. 



2. Identifiabilite 

2.1 . Equations differentielles 

Renvoyons a fSl pour des rappels sur les corps, differentiels ou non. 

Soit fco le corps de base de caracteristique nulle, Q par exemple. Soit fco(O) le corps 
engendre par un ensemble fini 8 = {6*1 , . . . , 6*^} de parametres inconnus. Soit k la cloture 
algebrique de A;o(0)- Introduisons le corps fc(s) des fractions rationnelles en I'indetermi- 
nee s, que Ton munit d'une structure de corps differentiel grace a la derivation ^ (les 



4. Lobry [31] a ecrit un pamphlet edifiant sur I'histoire « agitee » et la reception « houleuse » de cette 
analyse, en depit (a cause ?) de sa beaute et de sa puissance inconstestables. Les propres deboires 
de l'auteur lui ont prouve que ce temoignage n'est pas exagere I 

5. Voir, par exemple, 1 1 4] QH [T9] IS] |22| |23] et leurs bibliographies. Des questions classiques sur 
I'identification parametrique, les observateurs, le diagnostic et I'attenuation de perturbations y re- 
Qoivent des solutions d'une grande simplicite conceptuelle et faciles a mettre en oeuvre en temps reel. 
Mentionnons aussi la commande sans modele 1 15j, particulierement prometteuse. 

6. Voir (Te) pour une premiere ebauche. 
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elements de kg, de 8 et, done, de k, sont des constantes). Tout signal x, x ^ 0, est sup- 
pose satisfaire une equation differentielle lineaire homogene, a coefficients dans k(s), et 
done appartenir a une extension de Picard-Vessiot de k{s). 



Remarque 6 II suffit pour se convaincre de I'existence d'une telle equation homogene de 
de river les deux membres de (|7} suffisamment defois par rapport a s. 

L'anneau non commutatif k{s)[-j^] des operateurs differentiels lineaires 

finie 

est principal a droite et a gauche (cf. [361). Le ^(s)[^]-module a gauche engendre par x 
et 1 est un module de torsion (cf. Il36l ). et, done, un fc(s)-espace vectoriel de dimension 
finie, n + 1, n > 0. D'oii le resultat suivant qui semble nouveau (cf. If8l l4ll ) : 



Proposition 2.1 II existe un entier minimal n > 0, tel que x satisfait I 'equation differen- 
tielle lineaire, d'ordre n, non necessairement homogene, 

oii les polyndmes p, qq, . . . , Qn G k[s] sont premiers entre eux. Cette equation, dite mini- 
male, est unique a un coefficient multiplicatif constant non nul pres. 



2.2. Identifiabilite lineaire projective 

Rappelons que I'ensemble G = {di, . . . , 9g} de parametres est dit (cf. 112111231 ') 

- lineairement identifiable si, et selement si. 




ou 

- les entrees des matrices 21, carree g x g, et colonne g x 1, appartiennent a 

spanfc^(^)[_^](l,x); 

- det(2t) ^ 0. 

- projectivement lineairement identifiable si, et seulement si, 

- il existe un parametre, 9i par exemple, non nul, 

- I'ensemble {|^, . . . , est lineairement identifiable. 
Reecrivons ^ sous la forme suivante : 
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oil les + 1 coefficients a^,y et les M coefficients appartiennent a k. La matrice carree 
9Jl d'ordre N + M + 1, dont la ligne, < ^ < iV + M, est 

est singuliere d'apres dU et (|4|i. La minimalite de (|2|i permet de demontrer selon des 
techniques bien connues sur le rang du wronskien (cf. ^ |4T|) que le rang de 2Jl est 
N + M. II en decoule : 

Theoreme 2.2 Les coefficients a^i, et 6^ de d?} sont projectivement lineairement identi- 
fiables. 

Corollaire 2.3 Posons x = oii les polyndmes p,q (z k[s] sont premiers entre eux. 
Alors, les coefficients de p et q sont projectivement lineairement identifiables. 

II est loisible de supposer 1' ensemble des parametres inconnus 6 — {9i, . . . ,9g} stricte- 
ment inclus dans celui des coefficients a^^ et 6^ de (|4|, et done lineairement identifiable. 



3. Perturbations et estimateurs 

Avec une perturbation additive w le capteur fournit non pas x mais x + w. Soient 

- R = A:o(0)[s](/co[s])"^ I'anneau localise (cf. [29]) des fractions rationnelles a nu- 
merateurs dans fco(8)[s] et denominteurs dans fco[s]. 

- R[^] I'anneau non commutatif des operateurs differentiels lineaires a coefficients 
dans R. 

On obtient, a partir de la 

Proposition 3.1 Les parametres inconnus verifient 

( ei 

21 



= 05 + £ (6) 



V 



oil les entrees de €, matrice colonne g x 1, appartiennent a span^j_<i_] (w). 

On appelle (|6]l un estimateur. II est dit strictement polynomial en i si, et seulement si, 
toutes les fractions rationnelles en s, rencontrees dans les coefficients des matrices 21, 5B, 
€ de (|6]l, sont des polynomes en - sans termes constants. On peut toujours s'y ramener 
en multipliant les deux membres de ^ par une fraction rationnelle de kQ{s) convenable. 
On aboutit, alors, dans le domaine temporel, aux estimateurs consideres en [,1 IJ . si Ton 
suppose I'analyticite du signal : 



S{t) i[0,Ut) - 0,,) - C ' T^w{n)dndT2 ...dTk L=l,...,g 

Cm- "'0 "'O "'O 

(7) 

oil 
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- c est une constante, 

- [0, t] est lafenetre d' estimation, de largeur t, 

- 5{t) est une fonction analytique, appelee diviseur, nulle en 0, 

- [0]e(i) est I'estimee de 9 en t. 



4. Bruits 

Renvoyons a ll42l et ||9] pour la terminologie de F analyse non standard, deja utiUsee 
en flTl. Les propositions 14. ll et 14. 2l ci-dessous affinent la proposition 3.2 de |fTT|, oil les 
estimations sont obtenues en temps limite, « court » en pratique. 

4.1 . Sinusoides hautes frequences 

La perturbation du §[3]est de la forme 

M 

1=1 

ou 

- M est un entier limite standard, 

- les frequences 17, > sont des constantes ilHmitees, 

- les amplitudes sont des constantes, limitees ou non, 

- les phases (yS, , Q < < 27r, sont des constantes. 

Si les quotients ^ sont infinitesimaux, c'est un bruit centre, c'est-a-dire de moyenne 
nulle, au sens de Iflll . Des manipulations elementaires des integrates iterees O conduisent 
a la 

Proposition 4.1 Si 

- les quotients ^ sont infinitesimaux, et, en particulier, si les sont limites, 

- la largeur de la fenetre d' estimation est limitee et n ' appartient pas au halo d'un zero 
du diviseur, 

les estimees des parametres inconnus, obtenues grace a 0, appartiennent aux halos de 
leurs vraies valeurs. II n'en va plus de mime si I'un des quotients ^ est appreciable. 

Remarque 7 II existe des valeurs illimitees des amplitudes A^, v^HT par exemple, telles 
que les estimees precedentes appartiennent aux halos des vraies valeurs. 

4.2. Bruits blancs 

Designons par *N, *R les extensions non standard de N, M. Remplagons I'intervalle 
[0, 1] C M par I'ensemble hyperfini I {0, . . . , 1}, oil TV G *N est illimite. Un 
bruit blanc centre est une fonction w : I *M, l i-^ w{l) — An{L), oh 

- r amplitude A E *M est constante. 
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a2 

- le quotient ^ est limite, 

- les n{i) sont des variables aleatoires reelles, supposees centrees, de meme ecart-type 
1 normalise, et deux a deux independantes. 

Remarque 8 Cette definition non restreinte au cas gaussien, qu 'il convient de comparer d 
celle de fl^, precise [11 ] ; elle est inspiree de publications d'ingenieurs sur le bruit blanc 
en temps discret (voir, par exemple, [39]). Elle clarifie, a la maniere de [37], I'approche 
en temps continu usuelle dans les manuels de traitement du signal ( voir, a ce sujet, /E] 
ED \40]l et leurs bibliographies). Rappelons que cette approche continue est basee, en 
general, sur V analyse de Fourier et renvoyons, a ce sujet, a l ilOL Mentionnons, enfin, les 
travaux de 4251 1261/ . bases sur V analyse fonctionnelle. 

Remarque 9 Un pas supplementaire, inutile ici pour nos besoins, consisterait a rempla- 
cer, comme en H37]l . les variables aleatoires n{i) par des analogues « discrets ». 

Comme au § 14.11 il vient : 

Proposition 4.2 Si 

.2 

- le quotient — est infinitesimal, et, en particulier, si A est limite, 

- la largeur t, t £ I, de lafenetre d' estimation n'appartient pas au halo d'un zero du 
diviseur, 

les estimees des parametres inconnus, obtenues grace a appartiennent presque su- 
rement aux halos de leurs vraies valeurs. II n 'en va plus de meme si le quotient est 
appreciable. 

Remarque 10 II existe des valeurs illimitees de A, \/N par exemple, telles que les esti- 
mees precedentes appartiennent presque surement aux halos des vraies valeurs. 

Remarque 11 II est loisible de remplacer I'independance de et n{i,'), l ^ i! , par le 
fait que I'esperance du produit est infinitesimale. 
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